Strong Coupling Problem with Time- Varying Sound Speed 



Austin Joyce and Justin Khoury 

Center for Particle Cosmology, Department of Physics and Astronomy, 

,-H University of Pennsylvania, Philadelphia, PA 19104 

O 

Qh Abstract 



(N 



For a single scalar field with unit sound speed minimally coupled to Einstein gravity, there 

are exactly three distinct cosmological solutions which produce a scale invariant spectrum 

1-^ of curvature perturbations in a dynamical attractor background, assuming vacuum initial 

o, conditions: slow-roll inflation; a slowly contracting adiabatic ekpyrotic phase, described by 

(— I a rapidly-varying equation of state; and an adiabatic ekpyrotic phase on a slowly expanding 

background. Of these three, only inflation remains weakly coupled over a wide range of 

%^ modes, while the other scenarios can produce at most 12 e-folds of scale invariant and gaussian 

I ^^..^.„^^^^^^^^^.^...^.^^. 

lO in time affects this classification. While in the presence of a variable sound speed there 

CO 

are many more scenarios which are scale invariant at the level of the two-point function, 

^^ they generically suffer from strong coupling problems similar to those in the canonical case. 

,_H There is, however, an exceptional case with superluminal sound speed, which suppresses 

^ non-gaussianities and somewhat alleviates strong coupling issues. We focus on a particular 

^ realization of this limit and show these scenarios are constrained and only able to produce 

C^ at most 28 e-folds of scale invariant and gaussian perturbations. A similar bound should 

hold more generally — the condition results from the combined requirements of matching 
the observed amplitude of curvature perturbations, demanding that the Hubble parameter 
remain sub-Planckian and keeping non-gaussianities under control. We therefore conclude 
that inflation remains the unique cosmological scenario, assuming a single degree of freedom 
on an attractor background, capable of producing arbitrarily many scale invariant modes 
while remaining weakly coupled. Alternative mechanisms must inevitably be unstable or 
rely on multiple degrees of freedom. 
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1 Introduction 

There is clear evidence for a nearly Harrison-Zel'dovich and gaussian spectrum of perturbations 
from cosmic microwave background and large scale structure observations. One of the most 
pressing questions in early-universe cosmology is to understand the mechanism underlying the 
generation of this primordial spectrum. Inflation [1] is one explanation, but it is important to 
understand what other scenarios predict a spectrum of curvature perturbations consistent with 
observations. This question has motivated much research into alternative theories; for instance the 
pre-big bang scenario of [2] — [4], string gas cosmology [5] — [10] and ekpyrotic cosmology [11] — [43]. 
Production of a scale invariant spectrum of fluctuations by itself is not enough — it is also 
highly desirable for the background solution to be an attractor} Technically, this is achieved by 
demanding that the curvature perturbation on uniform-density hypersurfaces, C, goes to a constant 
in the long wavelength limit A; — )■ 0. In this limit, (^ ^ 5a/ a is interpreted as a constant perturbation 
of the scale factor, which may therefore be absorbed locally by a spatial diffeomorphism [46]. 

For models involving a single, canonical scalar field {i.e. with unit sound speed, Cg = 1) minimally 
coupled to Einstein gravity, it has been shown recently that there are only three independent 
cosmological solutions which produce a scale invariant spectrum of curvature perturbations on 
an attractor background [47, 48], assuming vacuum initial conditions. See [49] — [51] for related 
work. The most well-known of these solutions is of course inflation [1], which relies on exponential 
expansion of the background with e = —H/H"^ ~ 0. More recently, the adiabatic ekpyrotic [38, 39] 
scenario has been proposed, in which a scale invariant spectrum is produced by a rapidly evolving 
equation of state e ~ 1/t^ on a slowly contracting background. The third solution can be viewed 
as a variant of adiabatic ekpyrosis, where curvature perturbations are again sourced by a rapidly 
changing equation of state, but this time on a slowly expanding background [40]. At the level of 
the two-point function, these three scenarios yield indistinguishable power spectra. 

However, the degeneracy is broken at the three-point level. The non-inflationary solutions have 
strongly scale dependent non-gaussianities [39, 47], which can be traced to the rapid growth of 



""^The requirement that the background be an attractor may not be essential. Indeed, there are scenarios where 
instabihties play a crucial role and have important consequences, for example the matter bounce scenario [44] in 
the single-field case, as well as the curvaton mechanism and the phoenix universe [45] in the multi-field case. See 
[48] for a detailed discussion of single-field, non-attractor scenarios. 
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the equation of state parameter. In these models, /nl ~ k grows rapidly at small scales and 
perturbative control is lost when /nlC ~ 1- This difficulty can be avoided by suitably modifying 
the potential so that ( becomes much smaller on small scales. But this in turn restricts the range of 
scale invariant and gaussian modes to about 5 decades (~ 10^) or ~ 12 e-folds in fc-space [39, 47]. 
As a result, with Cg = 1 and attractor background, inflation is the unique single field mechanism 
capable of producing many decades of scale invariant and gaussian perturbations. (Of course, 
as a theory of the early universe inflation must still surmount some foundational issues, such as 
the measure problem and low-entropy initial conditions [52]. Here we leave aside these critical 
questions and note that inflation, viewed as a mechanism for generating density perturbations, 
remains weakly-coupled over a large range of modes.) 

In this paper we generalize the analysis to the case of time- varying sound speed, Cgit), as obtained, 
for instance, with non-canonical scalar fields. With Einstein gravity plus a single degree of freedom, 
the sound speed is the only remaining knob at our disposal.^ As shown in [53], allowing for Cs{t) 
greatly broadens the realm of allowed cosmologies that yield a scale invariant power spectrum. In 
particular, any cosmology with constant equation of state can be made scale invariant by suitably 
choosing the evolution of the sound speed. In this work we show that non-gaussianities impose 
stringent constraints on the allowed cosmologies. Our analysis is very general and applies to 
arbitrary time-dependent e(t) and Cs{t), with the only restriction that the null energy condition 
be satisfied: e > 0. 

We begin by reviewing how the time-dependence of the sound speed results in an effective cos- 
mological background for the curvature perturbation, as was ffist shown in [53]. In this effective 
background, which depends both on the evolution of the scale factor and the sound speed, ( 
propagates at the speed of light. We derive a consistency equation that the scale factor and the 
sound speed must satisfy in order to have scale invariance at the two-point level. In the spirit 
of [53], given an evolution for the scale factor, solving this equation gives a suitable evolution for 
the sound speed for which ( has a scale invariant two-point function on an attractor background. 
This shows that a time-dependent sound speed vastly increases the degeneracy at the two-point 
level. 



^One could also consider alternative theories of gravity. Our analysis applies to any theory of gravity which 
admits an Einstein frame description in terms of some field variables, such as generic scalar tensor theories. 
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As in the canonical case, this degeneracy is generically broken by the three-point function. In 
particular, if the three-point function is strongly scale-dependent, we generically expect the theory 
to become strongly coupled either in the infrared (IR) or in the ultraviolet (UV). To avoid such 
perturbative breakdown, we demand that certain contributions to the three-point function be 
scale invariant. This turns out to be extremely restrictive: we show that slow-roll inflation is the 
unique cosmology with this property. Conversely, if the three-point function is not scale invariant, 
then non-gaussianities will increase rapidly with scale, resulting in a finite range (< 10^ modes) 
of perturbations consistent with observations, as in the canonical case. This is a remarkable fact; 
it is extremely surprising, in light of the vast degeneracy afforded by a variable sound speed, that 
slow-roll inflation should be the unique possibility. 

There is, naively, a way to avoid our argument; since all the vertices in the cubic action are 
proportional to l/cj., it may be possible to have a strongly scale dependent three-point function 
while remaining weakly-coupled by taking c^ ^ 1. There has recently been much interesting work 
done related to this idea. A concrete example of this limit is given by the scenario considered 
in [55] — [60], based on earlier ideas in [61] — [64]. In these scenarios, the sound speed is assumed 
^ 1 in the very early universe and decreases in time. It should be noted that theories with 
superluminal propagation cannot be UV completed by a theory with an analytic S-matrix, for 
example local quantum field theory or perturbative string theory [54]. However, for the sake of 
generality, we remain agnostic about the nature of any UV completion of the theories we consider. 

We review the particular superluminal mechanism where e and e^ are constant and show that such 
cosmologies are also tightly constrained by strong coupling considerations. Again in these cases, 
only a finite range of scale invariant modes can be generated without hitting strong coupling. In 
order to understand this constraint, recall that the ( power spectrum is given by 

It is immediately clear that at early times, when Cs ^ 1, there is a tension between keeping 
H/Mp\ < 1 and maintaining the observed amplitude for density perturbations. For e ~ 0{1), this 
keeps us from taking Cg to be too large, which bounds the range of scale invariant modes to about 
8 decades (~ 10^). This tension may be alleviated by taking e to be exponentially large, but this 
causes the amplitude of non-gaussianities to grow. Because Cs is decreasing in time, demanding 
that the theory remain weakly coupled imposes a bound of about 12 decades (~ 10^^) on mode 
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production. Though this range is significantly larger than the ~ 10^ viable modes in the canonical 
case, we still have to account for the coincidence that these overlap with the observable window. 

Our results have some overlap with that of Baumann et al. [48], who reached similar conclusions. 
However, our analysis is much more general — in the variable sound speed case, Baumann et al 
considered only sub-luminal Cs and nearly constant e. In contrast, we allow arbitrary evolution for 
e and Cg — including the superluminal regime — in Section 3 and analyze in detail the constant 
e superluminal scenarios in Section 4. 



2 Scale Invariance with Variable Sound Speed 

We begin by considering how a varying sound speed affects the classification of scale invariance at 
the two-point level recently undertaken in [47, 48] . We make no assumptions about the underlying 
dynamics, only that they may be well modeled by a perfect fluid. Perturbing around a Friedmann- 
Robertson- Walker (FRW) background in co-moving gauge, where 5p = and hij = a(t)e^'^6ij, 
the quadratic action for ( is given by [65] 
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^2 = M^i / d^xdr z' 



f^] --im 



(2) 



where z = a^jcg^ and r denotes conformal time, adr = dt. This action is familiar from canonical 
single field models, except for the sound speed factor multiplying the spatial gradient term and 
appearing in the measure factor. In order to eliminate this complication, following [53] we define 
the sound horizon time coordinate by dy = c^dr. (Note that when Cg = constant, the variable y 
measures the size of the sound horizon.) Additionally, we define 

q = ^CgZ = -— . (3) 



(4) 



In terms of these new variables, the quadratic ( action takes the familiar form 



S2 = M^, f d^xdy q^ 



C - {dCY 



where primes indicate derivatives with respect to y. The virtue of this change of variables is 
manifest — ( now propagates luminally, but in effective cosmological background defined both by 
the scale factor and the sound speed. 
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The mode functions of the canonically normahzed field, v = \/2Mpi q ■ C, obey the famihar-looking 
equation of motion 

vl +(k'-^-^y, = 0. (5) 

This has exactly the same form as the corresponding equation in the canonical case, but with 

q replacing z, a manifestation of the curvature perturbation evolving in an effective geometry. 

Assuming the usual adiabatic vacuum, it is well known that (5) will yield a scale invariant spectrum 

of perturbations provided that 

o" 2 

- = 4 ■ (6) 

q y^ 

Note that modes freeze out when k\y\ ~ 1, which corresponds to sound-horizon crossing in the 

constant Cg case, hence we take — oo < y < 0. The solution for the mode functions is then 



which describes a scale invariant spectrum, Vk ~ /c"'^^^, in the limit y — ;■ 0. 

Equation (6) has two solutions, q ~ l/(— y) and q ~ y^, but only the former describes a background 
which is a dynamical attractor. To see this, note that in the long wavelength {k — > 0) limit we 
have the following expression for the power spectrum of the solution (7) 

n = Afc'iar~4^, (8) 

which is indeed independent of k. When q ~ l/(— y), (^ -^ constant outside the horizon, indicating 
perturbative stability [46]. The other solution q ~ y^, however, implies that C, grows outside the 
horizon, C, ~ y~^, signaling that the background is unstable. Since we are interested in attractor 
backgrounds, we henceforth ignore the q ^ y"^ solution. Recalling the definition of q, the condition 
for scale invariance in an attractor background may therefore be succinctly expressed as 

where /3 is an arbitrary (positive) constant. It is important to note that a and e are not independent 
degrees of freedom, but are related by e = —H/H^. Changing time variables to y, this relation 
becomes 

di c.ra^y ^^^^ 



dt H a \ c.^a' 
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Using the condition (9) for scale invariance, we can rewrite this to obtain the master equation 

4-)'-^.- (11) 

This equation ensures a scale invariant spectrum on an attractor background. As noted in [48], 
in the case where c^ = constant, this equation may be recast as a particular instance of the 
generalized Emden-Fowler equation. 

For completeness, we review the results of [47, 48] . In the case of constant sound speed (without 
loss of generality, we may take Cg = 1), there are three distinct scale invariant solutions: 

• Inflation is a solution where the scale factor grows as ai^f ~ l/(— r) and the equation of 
state parameter is constant e^f ^1 [1]. To check that this is in fact scale invariant, we note 
that gfjjf ~ 1/''"^! where y ~ r because Cg is constant. 

• Adiabatic Ekpyrosis is a solution where the equation of state parameter varies rapidly, 
eek ~ V''"^5 while the background remains nearly static, Ogk ~ 1- Again, we can check 
that this gives a scale invariant spectrum g^j^ ~ V'^^- ^^ fact, this corresponds to two dis- 
tinct solutions, one where the background is slowly contracting [38, 39] and one where the 
background is slowly expanding [40]. It is important to note that in these scenarios modes 
freeze out on sub-Hubble scales and are subsequently pushed outside the horizon during a 
contracting ekpyrotic phase with constant e ^ 1. 

Returning to the general case, given any evolution for a we can find an evolution of Cg that will 
make the spectrum of perturbations scale invariant by solving (11). Alternatively, specifying a 
relation between the evolution of Cg and a is sufficient to determine the evolution. As a result, we 
see that there is an enormous amount of degeneracy at the two-point level. 

An excellent illustration of this degeneracy is the case of cosmologies with constant e. With 
constant Cg, as reviewed above inflation is the only solution that has constant e. But for more 
general sound speed, there is a power-law evolution for Cg that yields a scale invariant spectrum 
for arbitrary positive values of e. Indeed, constancy of e and e^ = Cg/Hcg is sufficient to deduce 
the scaling solutions 

a ~ (— i/)'+'''-i , Cg ~ (—1/)"+"'"^ . (12) 



Inserting these expressions into (11), we find that the solution is scale invariant for e^ = — 2e, in 
agreement with [53]. 



3 The Cubic Action 



Non-gaussianities offer a powerful tool for differentiating between the different cosmologies with 
degenerate power spectra. Since the precise form of the cubic action depends the underlying 
physics, we must choose to parameterize the microphysics in some way. A convenient and quite 
general choice is to consider a non-canonical scalar field 0, described by a P{X, (p) lagrangian 



S 



j4 / 

d Xy/—g 



^^ + P(.Y. 



(13) 



where X = ~^g^'^df^(f)du(f). It is straightforward to show that the energy density and sound speed 
are related to the choice of P by [65] 

2 P,x 



p = 2XPx - P 



P.X + 2XP vx 



(14) 



The cubic action for ( is derived in [66, 67, 68]. Making the transformation to the sound horizon 
time variable dy = c^dr, the action takes the form, up to a field redefinition. 



S. 



d^xdy 



-act <^ S 1 



+ 



a^e 



c 
2e, + l 



|)+2A}g + f (e-3 + 3c^)CC'^ 



D C (dCY 



2^2 



2a'e 



ci 



■C'^C 



V2 



(15) 



a e ( 7] 
'2^ 



ec.§ec'4c.^'^--^''' 



-VV, 
4c? M V2 



where t] = H ^dlne/dt, and 



A = X^Pyx + ^X^P. 



XXX 



S = XP V + 2XP.X 



H\ 



(16) 



Although we focus on a particular class of microphysical models, namely P(X, 0) theories, the 
model-dependence of the action is encoded only in A. All other vertices in the cubic action are 
functions of the scale factor and the sound speed. For example, for a DBI action, the C,'^ term 



in (15) vanishes identically [53, 68]. Since the form of this first term will not be material to our 
arguments, our analysis even at the cubic level is rather general, but there may be some potential 
model-dependent effects from the ('^ vertex which we have not considered. 

To estimate non-gaussianities, a useful approximation is the horizon-crossing approximation, 
whereby /nl is estimated by 



/] 



NL 



C-A 



(17) 

k\y\=l 



Here £2 and £3 are terms in the quadratic and cubic lagrangians, respectively. Since temporal 
and spatial gradients are comparable at horizon crossing {dy ~ 9j ~ /c), we may trade them freely 
in (17). The horizon-crossing approximation generally offers a good estimate of /nl since modes are 
in their ground state at early times — when they are far inside the horizon — and become constant 
outside the horizon. We therefore expect non-gaussianities to peak around horizon crossing.'^ 

At a classical level, perturbations are highly non-gaussian for /nlC ^ 1; corresponding to C3/C2 ^ 1, 
and classical perturbation theory breaks down. At a quantum level, the right hand side of (17) 
also offers an estimate for the magnitude of loop corrections to the two-point function [69]. Thus, 
classical and quantum perturbation theory break down, and the theory becomes strongly coupled, 
whenever 

•-2 



cr' 



or /nlC ~ 1- This is the same strong coupling criterion used in [48]. 

In particular, if /nl is strongly scale dependent, then the growth of non-gaussianities will gener- 
ically lead to a breakdown of perturbation theory either in the IR or in the UV. (As mentioned 
earlier, an exceptional case is the limit c^ ^ 1, which will be treated separately in Section 4.) 
This expectation is borne out by the analysis of the canonical case [39, 40, 47]. To avoid strong 
coupling, therefore, we demand that /nl be approximately scale invariant. 



•^An important exception is the adiabatic ekpyrotic solution, where the e^ contributions peak at late times, 
well after horizon crossing [39]. Although ( goes to a constant outside the horizon, the rapid growth of the 
vertex, e^ ~ 1/^^, overwhelms the suppression from (^ derivatives becoming small. Thus, the horizon-crossing 
approximation is a conservative estimate of non-gaussianities. 
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Amongst the terms in the cubic action is the vertex^ 



^3 D / ^'xdy '^dC%C' . (19) 



2c3 "V2 

Evaluating this vertex at horizon-crossing, we find that its non-gaussian contribution is 



.,3 2c3 ^^ y: 



(20) 



k\y\=i 

Substituting the condition (9) for scale invariance at the two-point level, a?e/cs ~ 1/y^, this 
reduces to 

„g3 

/N] 



a4y4 



(21) 
fc|j/|=i 

Now, in order for the full three-point function to be scale invariant, a necessary condition is that 

the contribution from this vertex be scale invariant, barring miraculous cancellations. This implies 

that the scale factor must be growing as 

.^^, (22, 

which corresponds to an effective de Sitter geometry. Remarkably, simply demanding scale invari- 
ance of the two- and three-point correlation functions, without any consideration of the indepen- 
dent dynamics of a and c^, has led us to focus on backgrounds that are effectively de Sitter, albeit 
in terms of the y variable. Thus the question becomes — is it possible to have inflation without 
inflation? By this, we mean, is there an evolution where the modes see an effective de Sitter space 
in terms of the y variable but for which the true geometry is far from de Sitter? Unfortunately, 
the answer appears to be no, as we now argue. 

With a{y) ~ l/(— y), (9) immediately implies that e/cg = 7, where 7 is an arbitrary (positive) 
constant. This is all we need to solve (11), with the result 

-f log (y/y) log (y/y) 

where < |y| < \y\. Both e and Cg start out inflnite and decrease rapidly to zero. By construction, 
this solution is scale invariant at the two-point level and the aforementioned three-point vertex 
is also scale invariant. At flrst sight, we might expect that this solution is far from de Sitter 



''This vertex is the leading contribution to non-gaussianity in the adiabatic ekpyrotic scenarios [38, 39, 40, 47]. 
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because e ^ 1 initially. However, because e is decreasing so rapidly, by the time \y\ < e^^||/|; e is 
already less than unity, indicating an inflationary spacetime. As such, this solution is only a small 
deformation away from the de Sitter geometry, specifically only about one e-fold of evolution is 
non-inflationary. 

It is worth noting that while the particular vertex (19) we have considered yields a time- independent 
(and therefore scale invariant) non-gaussian contribution once (22) is imposed, some of the other 
vertices vary logarithmically with y, such as the C,'dC,^ vertex. But this slow growth as a func- 
tion of scale is of course acceptable from the viewpoint of avoiding large non-gaussianities and 
the breakdown of perturbation theory. 

The analysis thus far strongly suggests that if we want scale invariance at the level of the three- 
point function, we are forced to consider an evolution that is very close to de Sitter — that is, 
the cosmological background must inflate. The intuition for this statement is that every term in 
the cubic action is multiplied by powers of e. In the case of slow-roll inflation, we take e ^ 1 
and therefore get a nearly gaussian spectrum. Note that there is another factor common to every 
term, namely 1/c^. One might therefore suspect that a nearly gaussian spectrum could also result 
from taking c^ :» 1, even with a strongly scale dependent /nl- This will be treated separately in 
Section 4, where we will show that other considerations should greatly constrain this case as well. 

One may wonder whether the assumption of scale invariant three-point function is not too strin- 
gent. A possible loophole would be to have non-gaussianities peak in the IR, where the onset 
of strong coupling would lie on unobservably large scales. However it is clear from the above 
analysis that this would require a scale factor that expands more rapidly than a ~ l/y, which is 
super-inflationary in terms of the y variable. We have studied this case by generalizing the above 
analysis to power-law dependencies for /nl, and found that inflation (in terms of cosmic time) 
remains the only solution. 



4 Superluminal Cs Cosmologies 

Given the above arguments, it seems as though it is impossible to have small non-gaussianities over 
a broad range of modes without invoking a period of inflation. However, there is an illuminating 
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exceptional case which we now treat separately. As we noted in passing earlier, every term in the 
cubic action has a factor of l/c^ multiplying it. In analogy with slow-roll inflation, we may imagine 
controlling the size of the three-point function by making this factor very small, corresponding to 
superluminal Cg. In principle, this should allow for the production of modes with an extremely 
gaussian spectrum over a broad range in /c-space [58]. In practice, however, we will find that only 
a finite range of scale invariant and gaussian modes can be generated by such mechanisms. 

Superluminal scenarios have been proposed recently [55] — [60], building on previous work on de- 
caying sound speed cosmologies [61] — [64]. These models take advantage of having an extremely 
large speed of sound at early times which suppresses non-gaussianities on large scales. By virtue 
of the fact that the sound horizon at early times is much larger than the co-moving Hubble radius, 
the horizon problem is addressed. The large sound speed decays as time goes on, causing the co- 
moving sound horizon to shrink in time. See Fig. 1. In the same manner that a shrinking Hubble 
horizon causes modes to freeze out in inflation, modes freeze out when they exit the sound hori- 
zon in this case. Additionally, the flatness and horizon problems are effectively decoupled in this 
scenario. Of course, these models manifestly rely on superluminality, which, as mentioned before, 
makes them unsuitable for UV completion within a theory that is local in the usual sense [54] . 

Current models in the literature [55, 56, 60] are cast in terms of constant e cosmologies mentioned 
earlier and studied extensively in [53]. In these scenarios, we assume scaling solutions for a and 
Cs given by (12): 



^■"■"'; o,= UV"-. (24) 

yij \yiJ 

where j/i is the time when Cg = 1. Using the freedom to rescale spatial coordinates, we have set 
a = 1 at this time also. The Hubble parameter is given by 



H=^-^ Tt(-) • (25) 

yi(e + e,-l) V^i/ 



The contribution to /nl given by (20) scales as 



rsj 



eV 1 



^^^'^v^; '^7 



r\j 



k' , (26) 



k\y\=l 

which thus peaks on small scales. Despite this strong scale dependence, the range of modes 
accessible to observations can be highly gaussian by having c^ 3> 1, which amounts to making the 
coefficient in (26) very small. 
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Figure 1: The initial setup with the sound horizon far outside the Hubble horizon. As the y 
variable evolves, the comoving sound horizon decreases and the comoving Hubble horizon grows 
until the two meet at the coincidence scale denoted by the dashed line. 



First note that the Hubble horizon at mode freeze-out and the amplitude of the power spectrum 
have a simple relationship. Indeed, the power spectrum in the constant e case is given by [53] 



{e + es- 11 



H^ 






^vr^e 



CsMl, 



1 



(%i 



2e + ta 



(27) 



k\y\=i S^r^e yl 

which confirms the statement made earlier that exact scale invariance requires e^ = — 2e. We will 
henceforth assume that e^ = — 2e. It then follows from (27) that the Hubble parameter and the 
power spectrum are related by 



H' 



M' 



^Tlh 



1 + e 



^CsPc , 



(2^ 



where everything is evaluated at freeze-out. Staring at this expression, it is immediately clear that 
in the limit c^ ^ 1 a tension arises between matching the observed normalization of P,^ ~ 10^^° 
and keeping H < Mpi. 

To be completely general, suppose that the phase of decaying sound speed lasts from some initial 
time, yinitiah to some final time ygnai- Without loss of generality, we will assume that ysnai < Hi 
— there is no gain in having the decaying Cg phase lasts any longer, since the modes generated 
when y > yi will never be outside the Hubble horizon. Thus the range A^ of scale invariant modes 
is given by ^/initial = A^l/finai- Substituting this into (25), we find that the ratio of the Hubble 
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parameter at the beginning and end of the decaying c^ phase satisfies 

-'^initial / yinitial \ 



-f^final V 2/final / 

Meanwhile, from (28) we have 



N^^ . (29) 



10 °-— VC^Ifinal , (30) 



Mpi 1 + e 



where we have used the observational constraint Pr ~ 10 ^°. We therefore obtain 



^intial , ^-5 V87r2"e 



For consistency of classical gravity, we demand that the Hubble parameter be sub-Planckian at 
the onset the decaying sound speed phase, -^initial ^ ^pi- This gives the following upper bound 

on A^: 

Remarkably, this bound on the range of scale invariant modes follows merely by requiring the 
correct amplitude for perturbations and that the physics be sub-Planckian. So far we have said 
nothing about strong coupling issues. It turns out that the bound can be made more precise by 
considering non-gaussianities. 

As before, consider the contribution to /nl given by (26). Demanding that this be consistent with 
perturbation theory requires 

/^k-f-^1 ■10-^<l, (33) 



5 

which implies the lower bound c^lmt ^ 10^2 ■ e. Inserting this into (32), the bound on A^ becomes 



A^< 



25 



10^ /1 + e 



Vs^ 



(34) 



Of course we are interested in the case e > 1, for otherwise the universe would be infiating. The 
bound (34) is weakest for e ~ 1, which is sensible, since in this regime c^ has to be just barely 
superluminal in order to suppress non-gaussianities. In this case, (34) enforces iV ~ 10^^ ~ e^^ 
for the range of scale invariant modes in fc-space. While this is sufficient to account for microwave 
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background and large scale structure observations, it does represent a relatively narrow range 
within which the observable window must lie. 

To summarize, it appears that strong coupling considerations also greatly constrain scenarios based 
on superluminality. The bound on the allowed range of scale invariant modes is weaker than in the 
canonical case (A^ < 10^^ versus < 10^, respectively), but is nevertheless restrictive. The scenario 
must somehow explain the coincidence of why the observational window happens to fall within the 
limited scale invariant range. Although here we have focused on the particular realization (24), 
we expect that similar arguments will bound mode production in other realizations of the limit 
Cs >1. 



5 Concluding Remarks 

A fundamental objective of early universe cosmology is to explain the origin of the observed 
spectrum of density perturbations. The primordial perturbations inferred from observations are 
nearly scale invariant and highly gaussian. Any mechanism that purports to explain their origin 
should be capable of generating perturbations with the desired properties over many decades in 
fc-space. 

The present paper concludes a body of work, initiated in [47, 48], aimed at identifying all pos- 
sible cosmologies capable of generating scale invariant and gaussian perturbations. The analysis 
assumes a single degree of freedom (purely adiabatic perturbations), and a background which rep- 
resents a dynamical attractor. In the case of c^ = 1, it was shown in [47, 48] that only inflation can 
generate suitable perturbations over many decades in fc-space. Two other phases, corresponding 
to adiabatic ekpyrosis, can produce at most 10^ modes (or 12 e-folds) of scale invariant and gaus- 
sian modes, before hitting strong coupling. Baumann et al. [48] also considered a time-dependent 
sound speed, but their analysis was restricted to sub-luminal Cg and nearly constant e. 

In this paper, we have generalized the analysis by allowing for arbitrary time-dependent e(t) and 
Cs(t), with the only restriction that the null energy condition be satisfied: e > 0. At the two-point 
level, any cosmological background can in principle yield a scale invariant spectrum, provided that 
the sound speed evolves in a suitable way. To avoid strong coupling we have argued that the three- 
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point function should generically also be nearly scale invariant. By demanding that a particular 
contribution to /nl be scale invariant, we have found that inflation remains the unique solution. 
Consequently, non-infiationary solutions have strongly scale-dependent non-gaussianities, and 
therefore the range of modes consistent with observations is subject to a bound similar to the 
canonical case. 

One exceptional case is the limit of superluminal sound speed, which has attracted considerable 
attention recently [55] — [60]. For concreteness, we have focused on superluminal scenarios where 
both e and eg are constant in time. In this case, non-gaussianities are strongly scale dependent, but 
strong coupling is naively avoided by making Cs ^ 1. While this does alleviate the constraints, we 
have shown that the combined requirements of sub-Planckian Hubble parameter and controlled 
non-gaussianities imply a bound of at most 10^^ modes (or 28 e-folds) with nearly scale invariant 
and gaussian statistics. This bound is weaker than in the canonical case, but is nevertheless quite 
restrictive. Although we have not considered it in detail here, we expect arguments similar to 
those in Section 4 to bound mode production in more general superluminal scenarios. 

We conclude that, under the assumptions of single degree freedom, attractor background and 
adiabatic vacuum initial conditions, inflation is the only mechanism capable of generating suitable 
perturbations over a wide range of scales. Alternative mechanisms must inevitably rely either on 
an instability of the background, non-vacuum initial state, and/or multiple degrees of freedom. 

Note added: As this paper was being completed. Reference [70] appeared, which reaches conclu- 
sions similar to ours using different arguments. Specifically, the authors find that that in order 
to produce more than 7 e-folds of scale invariant modes in a monotonically expanding universe 
without inflation — barring super-Planckian energy densities — a superluminal speed of sound 
must be invoked. However, this result relies on the assumption of an expanding background, while 
our arguments do not require this constraint. 
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